Coupling of bonding and antibonding electron orbitals in double quantum dots by 

spin-orbit interaction 
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We perform a systematic exact diagonalization study of spin-orbit coupling effects for stationary 
few-electron states confined in quasi two-dimensional double quantum dots. We describe the spin- 
orbit-interaction induced coupling between bonding and antibonding orbitals and its consequences 
for magneto-optical absorption spectrum. The spin-orbit coupling for odd electron numbers (one, 
three) opens avoided crossings between low energy excited levels of opposite spin orientation and 
opposite spatial parity. For two-electrons the spin-orbit coupling allows for low-energy optical tran- 
sitions that are otherwise forbidden by spin and parity selection rules. We demonstrate that the 
energies of optical transitions can be significantly increased by an in-plane electric field but only for 
odd electron numbers. Occupation of single-electron orbitals and effects of spin-orbit coupling on 
electron distribution between the dots are also discussed. 

PACS numbers: 73. 21. La 



I. INTRODUCTION 

In a pair of quantum dot&ir— defined in semiconduct- 
ing medium the charge carriers form extended wave func- 
tions when their tunneling through the interdot barrier 
becomes effective enough. In vertically stacked quan- 
tum dots the extended electron and hole orbitals are 
probed by photoluminescence experiments in external 
electric fieldj^ The electron single-dot orbitals hybridize 
to bonding ground states similar to the ones found in nat- 
ural covalent molecules. Recent studies^ indicated that 
the hole in artificial molecules of self-assembled quantum 
dots behaves in a different manner forming an antibond- 
ing ground-state orbital. This peculiar behavior results^ 
of the spin-orbit (SO) coupling induced mixing of light 
and heavy hole states. 

In the present paper we study the mixing of bonding 
and antibonding electron orbitals that is induced by SO 
interaction in planar systems of laterally coupled quan- 
tum dots. The coupling between spatial and spin elec- 
tron degrees of freedom results from inversion asymmetry 
of the structure^ and/or the crystal lattice.— This asym- 
metry enters into the two-dimensional SO Hamiltonian 
which does not conserve the spatial parity and couples 
the electron spin-up bonding orbitals with spin-down an- 
tibonding orbitals. In order to indicate experimentally 
accessible consequences of this coupling we consider op- 
tical absorption spectra in the external magnetic field for 
up to three confined electrons. In parabolic quantum dots 
the spin-orbit coupling introduces a distinct dependence 
of the far infrared magneto-optical absorption spectra on 
the number of confined clcctronSf^ We find that the SO 
induced modification to the absorption spectra of double 
dots are qualitatively different for even and odd electron 
numbers. 

Laterally coupled quantum dots^iiS are considered can- 
didates for realization of a quantum gate working on elec- 



tron spins^ since the height / width of the interdot barrier 
can be tuned by external voltages which is essential for 
the control of the spin exchange between the electrons 
confined in adjacent dots. The idea of the spin exchange 
motivated a number of theoretical investigations on the 
properties of electron systems in laterally coupled quan- 
tum dotsjii"— 



The SO interaction is one of the issues that are investi- 
gated in the context of spin-based quantum information 
processing.—"—!^"— The SO coupling allows for spin ma- 
nipulation by the spatial electron motion.—"— Moreover, 
it leads to the spin relaxation^SriSi mediated by phonons, 
leading to information decay and decoherence. Singlet- 
triplet induced avoided crossing of two-electron energy 
levels were observed in electron-transport spectroscopies 
for gated InAs nanowire quantum dots^^ as well as for 
double InAs quantum dotsi^ The exchange interaction 
between electrons confined in separate dots was found to 
contain an anisotropic component originating from the 
SO coupling^ which initially motivated a quest for spin 
processing procedure a^^'^^ minimizing its effects. Later 
on, proposals of using the asymmetry of the exchange 
interaction for construction of universal quantum gates 
that could work without single spin operations22i2i were 
formulated. Recently, it was demonstrated by the exact 
diagonalization thalfi^ the anisotropy of the exchange in- 
teraction - previously discussed within approximate ap- 
proaches - is in fact absent in zero magnetic field. There- 
fore, the treatment of spin-orbit coupling effects for dou- 
ble dots requires an exact diagonalization approach which 
we employ below. The SO coupled double quantum dots 
were so far studied by the exact diagonalization in Ref. 
[3^, which provides a detailed analysis of single-electron 
states and in Refs. [l8lf37| which deal with the electron 
pair in the context of the exchange interaction. 
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II. THEORY 



and the cubic 



We consider an effective mass singie-eiectron Haniilto- 
nian of tlic form: 



W{v] 



\2m* 

\g^BBaz + HsiA + Hbia, 



(1) 



where p = ?ik = —iKSI + eA, 1 is the identity ma- 
trix, W{y) stands for the potential, Hsia and Hbia 
introduce Rashbai^ (structure inversion asymmetry) and 
Dresselhaus^ (bulk inversion asymmetry) spin-orbit in- 
teractions. The vector potential is taken in the symmet- 
ric gauge A = -j(— 0). The Rashba and Dresselhaus 
SO interactions have the form 



and 



Hsia = aVW • (ct x k), 



Hbia = 7 [(Jxhooik^ - k^) + ayky{kl - kl) 



(2) 



(3) 



respectively. In Eqs. (2) and (3) a and 7 are bulk SO 
coupling constants, a's are Pauli matrices and x, y, z axes 
are oriented parallel to [100], [010] and [001] (growth) 
crystal directions, respectively. 

We assume that the confinement potential forming the 
quantum dot is separable into an in- plane Vc{x,y) and a 
growth direction Vz{z) components so that the potential 
appearing in the Hamiltonian (1) is 



Wir) = Vcix,y) + V4z) + \e\F 



(4) 



where F is the electric field vector (below we always take 
Fy = 0). In the following we adopt a two-dimensional ap- 
proximation assuming that the electrons occupy a frozen 
lowest-energy state of quantization in the growth direc- 
tion. The two-dimensional SO terms are obtained by 
averaging Hsia and Hbia over the wave function de- 
scribing the electron localization in the growth direction. 
The two-dimensional Rashba terms are usually^ sepa- 
rated into a diagonal 



H 



diag 
SIA 





'dV 




'dV , ,■ 




auz ^ 


dy 




k^ 



(5) 



and linear 



^SIA 



dW 
a{-^){crxky 
oz 



^ykx\ 



parts. In this formula the average gradient of the po- 
tential calculated for the wave function in the growth 
direction can be attributed to an effective z component 
of the electric field K. = rr^^r^). The two-dimensional 



Dresselhaus SO interaction contains the linear 



Tjlin 
^BIA 



(Tyky] 



(7) 



H^mA 7 [c^ykykl 



'^xkxky\ 



(8) 



terms. We assume that the quantum dot is made of 
Ino.5Gao.5As alloy for which we adopt the SO coupling 
constants a = 0.572 nm^ (after Ref. SO) and 7 = 32.2 
meVnm^ (after Ref. Issh . The other material parameters 
are taken as arithmetic average^ of InAs and GaAs, i.e. 
we use the electron effective mass m* = 0.0465tooj Lande 
factor g = —8.97 and dielectric constant e = 13.55. The 
considered large value of the g factor is in the order of 
the one found for in experimental sample o^^'^^ in which 
the SO coupling effects were studied. 

For the electron wave function in the growth direction 
identified with the ground-state of an infinite rectangular 
potential well of height d one obtains the two-dimensional 
linear Dresselhaus constant 7^^ = 7{fc^) = [see Eq. 
(7)] . In the bulk of our calculations we assume a minimal 
but still realistic value of d = 5.42 nm, for which 7^^ = 
10.8 meVnm. 

The in-plane confinement potential is taken in form 



Vc{x,y) = - 



Vo 



Vb 



(^ -[«])(' + [*]) 



(9) 



where Vb = 50 meV is the depth of the dots and Vh is 
the height of the interdot barrier. Wc assume ^ = 10 for 
which the potential profile has a form of a nearly rectan- 
gular potential well, where 2Rx = 90 nm and 2Ry ~ 40 
nm determine the size of the double dot in x and y di- 
rections respectively and 2i?f, = 10 nm is the thickness of 
the interdot barrier. Wc consider two values of the bar- 
rier height Vb = 10 mcV - for the double dot potential 
and Vfc = - for a single elongated dot. The elongated 
dot potential on the one hand corresponds to the limit 
case of strong interdot tunnel coupling and on the other 
it is close in geometry to the nanowire quantum dots, 
in which the spin-orbit coupling induced singlet-triplet 
avoided crossing was observed in a single-electron charg- 
ing cxperimenti^ The potential Vc is displayed in Fig. [T] 
for both the single and double dot. 

The single-electron eigenfunctions arc found by diago- 
nalization of the two-dimensional version of Hamiltonian 
(1) in a basis of multicenter Gaussian functions^ with 
embedded gauge invariance 



(6) i'n^Y^ c^,X. exp 



(r - Rk)2 ieB 



2a2 



+ ^{xYk^yXk) 



(10) 

where summation over k runs over centers of Gaussian 
R-k = (Afe, Yfc), s =t, i and Xs are eigenstates of Pauli 
matrix. The centers Rk are distributed on a rectangular 
mesh of 25 X 11 points spaced by Ax = Aj/ = 5.2 nm. The 
variationally optimal basis function parameter a = 4.7 
nm is used in the calculations. 
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FIG. 1: The shades of blue show the in-plane potential of a 
single dot (Vj, = - left column of plots) and of a double 
dot (14 — 10 meV - central and right columns). In the right 
column of plots an in-plane electric field of Fx = 0.5 kV/cm is 
included. Inside the light (darker) blue area the potential falls 
below -5 meV (-45 meV). The contours indicate the charge 
density for a single (top row) , two electrons (middle row) and 
three electrons (lowest row of plots) for B = 0. 



The eigenproblem of of A^-clcctron Hamiltonian 



N 



N 



H 



47rener 



(11) 



is solved using the configuration-interaction approach 
with a basis constructed of Slater determinants built of 
single-electron eigenfunctions (10) of SO-coupled Hamil- 
tonian. Convergence of the energies with a precision bet- 
ter than 1/icV is usually reached for inclusion of thirty 
one-electron eigenstates. 

The confinement potential (9) is symmetric with re- 
spect to the origin. In the present work the asymmetry 
effects are introduced by the in-plane electric field Fx- 
For Fx ^ and without SO coupling the stationary states 
possess a definite spatial parity with respect to point in- 
version P^/'"(— r) = ±-i/'"(r), where P is the inversion 
operator. The eigenvalue -1-1 corresponds to even parity 
states and the eigenvalue —1 to the odd parity states. 
When SO is introduced the spatial parity eigenvalue is 
no longer a good quantum number even for Fx = 0. For 
symmetric systems the SO coupled Hamiltonians com- 
mute with the operator Paz , which implies that the spin- 
up and spin-down components still possess definite but 
opposite spatial parities. We refer to Pa^ as the s-parity 
operator. Eigenstates of this operator with eigenvalue 
-1-1 (-1) are referred to as even (odd) s-parity states or 
for brevity s-even (s-odd) states. The even s-parity states 
have even-parity spin-up component and odd-parity spin- 
down component. 

We evaluate the optical absorption spectrum using the 
energies of stationary states and transition probabilities 
from state k to I that is proportional to the square of the 



N 

(*,|^(x,±zy,)l*;), 
i=i 



(12) 



where ^'fe is the A^-electron wave function for k-th Hamil- 
tonian (11) eigenstate and the signs ± correspond to op- 
posite circular polarization of the exciting light. The op- 
tical transitions conserve the electron spin and invert the 
spatial parity when it is a well-defined quantity. When 
the SO coupling is introduced the optical transitions can 
only occur between states of opposite s-parity. 




B[T] 

FIG. 2: Lowest single-electron energy- levels in function of the 
magnetic field for a single elongated dot (14 = 0) (a) and for 
a double dot {Vt = 10 meV) (b). Solid (dashed) lines corre- 
spond to the even (odd) s-parity. Black lines show the results 
without SO coupling. The blue curves show the results ob- 
tained with SO coupling for Fz = and the red curves in (a) 
for Fz = 188 kV/cm. The spin direction for the energy levels 
without SO coupling are marked with arrows. The thin verti- 
cal lines indicate allowed optical transitions from the ground- 
state. 
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FIG. 4: Contribution of even-parity orbitals (a) and average 
value of the z component of the spin (b) (in h units) in the 
lowest energy s-even and s-odd parity eigenstates. 



FIG. 3: Dashed curves show the potential confinement pro- 
file for Vb = 10 meV calculated for y = 0. (a) The spin-up 
components of the even s-parity ground state (red lines) and 
s-odd parity excited state (blue lines) in the absence of SO 
coupling (spin-down components exactly vanish), (b) Same 
as (a) but with SO coupling present. Spin-down components 
are presented in (c) The scale for the wave function is the 
same on all the plots, but in (c) the wave functions were mul- 
tiplied by 4. At right (left) panels an electric field is Fx — 0.5 
kV/cm (zero). 



III. RESULTS 
A. Single electron 

The single-electron spectrum for a single elongated dot 
and for the double dot is presented in Fig. [2] For B = 
the ground state and the first excited state are Kramers 
doublets. In each doublet we find one state of the odd 
s-parity and the other of the even s-parity. At i? = the 
electron in the ground-state (first-excite-state) doublet 
occupies predominantly a bonding (antibonding) orbital. 
With the solid (dashed) lines we plotted the even (odd) s- 
parity energy levels. Black lines show the results without 
SO coupling. The blue lines correspond to the case of 
SO coupling without the linear Rashba term {Hg}\), i.e. 
for Fz = 0. The red curves in Fig. ^a.) correspond to 
Fz = 188.8 kV/cm, for which the linear two-dimensional 
Rashba constant is as large as the linear two-dimensional 
Dresselhaus one. Beyond increased width of the avoided 
crossing no qualitative difference in the results obtained 
for these two values of F^ is found. Therefore, below we 
assume F^ = unless stated otherwise. 

For illustration of the double-dot wave functions we as- 
sumed a presence of a residual magnetic field B ~ 10/iT 
which lifts the doublet degeneracy and we chose the 
states of the ground and excited doublets that correspond 
to (sz) > 0. With the blue lines in Fig. [3] we plotted 
the spinor components of the even s-parity ground state 
which is bonding in its spin-up component with or with- 



out SO coupling. Its antibonding spin-down component 
appears when the SO coupling is introduced [Fig. EIc)]. 
The red lines in Fig. [3] correspond to the odd s-parity 
state of the excited doublet which is antibonding in the 
spin-up component. The SO coupling adds to this state 
a bonding spin-down component. 

In Fig. [5] one observes an avoided crossing of two ex- 
cited energy levels of the odd s-parity stemming of both 
the ground and the exited doublets. Without the SO 
coupling the energy level that goes up in the energy 
with growing magnetic field corresponds to the spin-down 
bonding orbital, and the one that goes down - to the spin- 
up antibonding orbital. The avoided crossing opened by 
the SO interaction is accompanied by spin and spatial 
parity mixing. 

For the single electron in ideally symmetric pair of 
dots (Fx ~ 0) there is a direct correspondence between 
the SO-coupling-induced mixing of both the spin states 
and the occupation of molecular orbitals of opposite spa- 
tial parity. The occupation of the even parity orbitals 
[oc(e)] is calculated as the norm of this component of the 
spinor that corresponds to the even parity state. Then 
the average value of the z-component of the electron 
spin is {sz) = ?i (oc(e) — i) for the even s-parity and 
[sz) = ^ (5 ~ oc(e)) for the odd s-parity states. Occupa- 
tion of the even parity orbitals and (s^) is for the double 
dot displayed in Fig. |4]as function of the magnetic field. 
The ground state at higher field becomes a pure bonding 
spin-up orbital. We notice that the values corresponding 
to the two odd s-parity energy levels interchange near 2 T 
which is related to the energy level anticrossing presented 
in Fig. [DJb). At the center of the avoided crossing these 
two energy levels correspond to (sz) = and bonding 
and antibonding orbitals are equally occupied. 

The discussed anticrossing of the odd s-parity energy 
levels leaves a clear signature on the optical absorption 
spectrum. The energy and probability of excitation from 
the ground-state are displayed in Fig. [SJ The ground- 
state has the even s-parity hence the absorption is only 
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FIG. 5: The dots show the low-energy absorption from the 
ground-state at = for the single dot (a) and for the 
coupled dots (b) (for the energy spectra see Fig. [2|. The area 
of the dots is proportional to the absorption probability. The 
black dots show the results without SO coupling. The full 
blue dots correspond to SO coupling with = 0. The open 
blue circles in (b) correspond to height of the dot d increased 
from 5.42 to 7.67 nm which amounts in a two-fold reduction of 
the 2D Dresselhaus constant. The red dots in (a) correspond 
to a strong linear Rashba coupling present = 188 kV/cm. 
The dashed grey line indicates the Zeeman splitting gfibB. 



allowed to the odd s-parity final state. The ground-state 
is nearly spin-up polarized (Fig. 21) and since electron 
spin is left unchanged during an optical transition the 
absorption goes to the s-odd state with spin-up orien- 
tation. When the avoided crossing is opened between 
the s-odd energy levels both of them possess a non-zero 
spin-up component and the optical transitions to both 
of them from the ground-state are allowed. Outside the 
avoided crossing the absorption spectra with or without 
SO coupling are similar. 

The energy range in which the SO-induced avoided 
crossing is observed in the absorption spectrum corre- 
sponds to far-infrared or microwave radiation in which 
cyclotron resonance experiments are performed.— One 
can increase the energy of the avoided crossing twice by 
applying an electric field of 0.5 kV/cm - see Fig. [BJb). 
In the presence of the electric field the electron in the 



ground-state is pushed to the left dot by Fj; > while 
the final state in the absorption process is mainly local- 
ized in the right dot [see Fig. [3{^b)]. The opposite shifts 
of the electron wave function in the initial and final states 
are translated by the electric field into an increased tran- 
sition energy [see Fig. |6l[a) for the energy splitting]. The 
obtained energy increase is accompanied by reduction of 
the SO-induced avoided crossing. 

Fig. El^b) shows also that for non-zero F the absorp- 
tion probabilities vanish at higher B. The separation of 
the initial and final states [Fig. Eljb)] by the electric field 
is enhanced when the magnetic field is applied, since the 
latter increases the localization of wave functions near 
the centers of the dots lifting the interdot tunnel cou- 
pling. In consequence - the ground-state becomes totally 
localized in one dot and the final state of the transition in 
the other. Vanishing overlap between the initial and final 
state wave function implies vanishing transition probabil- 
ity as calculated by formula (12). 



B. Electron pair 

In the absence of the magnetic field and without SO 
coupling the first excited state of the electron pair is spin 
triplet. For S = we find that the first excited state 
is threefold degenerate also with SO coupling present. 
This applies to both the single elongated dot [Fig. El^a)] 
and the double dot [Fig. [7Jb)]. Without SO coupling 
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FIG. 7: Two-electron energy spectrum for a single elongated 
dot Vt = (a) and for a couple of dots separated by Vt, = 10 
meV barrier (b). Black (blue) lines show the results without 
(with) SO coupling. For the results without SO coupling we 
added labels 5* for the singlet and T for the triplets (subscript 
denotes the sign of the z-component of the total spin) . Dashed 
(solid) curves correspond to odd (even) s-parity. Results were 
obtained for _Fa; = 0. 



the magnetic field induces a singlet-triplet ground-state 
transition near 1 T for the single dot and near 0.4 T 
for the double dot. The crossing singlet and triplet en- 
ergy levels have the same odd s-parity and an avoided 
crossing is opened between them when SO coupling is 
introduced. The calculated width of the avoided cross- 
ing is 0.18 and 0.07 meV for the single and double dot, 
respectively which is within the order of the ones found 
in experiments: 0.25 meV and 0.2 meV for the nanowire 
quantum dol— and for the double dots^^. 

For a symmetric system {F^ = 0) the optical transition 
from the ground-state can only go to the even s-parity 
eigenstate. In the absence of the spin-orbit coupling in 
the considered energy range only the triplet with zero z- 
component of the spin (Tq) has the required spatial parity 
to absorb photons. However, this absorption is excluded 
anyway on both sides of the singlet-triplet ground-state 
transition. For B below this transition the matrix el- 
ement (12) vanishes due to opposite symmetry of the 
spatial initial and final wave functions with respect to 
the electron interchange. For B above the singlet-triplet 
transition the ground-state (triplet with ~ h denoted 
as T+) and Tq states have the same symmetry with re- 
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FIG. 8: Two-electron ground-state absorption spectrum as a 
function of the magnetic field for SO coupled single dot (a,b) 
and double dot (c,d). Panels (a,c) correspond to Fx — and 
(b,d) to Fa; = 1 kV/cm. The area of the dots is proportional to 
the absorption probability. Transitions are denoted by labels 
of two-electron spin eigenstates which are found without SO 
coupling. Without SO coupling all the transitions presented 
in this figure are forbidden. 
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FIG. 9: Electron charge localized in the left dot as a function 
of the electric field for the double dot at B = 0. Results with 
and without SO coupling are not distinguishable. 



spect to the electron interchange, but the z-componcnts 
of the spin are different. Optical transitions between the 
states corresponding to energy levels presented in Fig. [7] 
are only allowed by the SO coupling. The calculated ab- 
sorption spectrum is shown in Fig. El For Fx = [Fig. 
[5{a,c)] the absorption probability grows with the mag- 
netic field after the singlet-triplet ground state avoided 
crossing. Then, the transition corresponds to — >■ Tq 
excitation in terms of states without SO coupling. When 
the electric field F^ is switched on [Fig. [Sl[b,d)] the parity 
selection rules no longer apply and we notice appearance 
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FIG. 10: Two-electron energy spectrum of the double dot is shown in (a,c,e) as a function of the electric field. Plots (b,d,f) 
indicate the charge localized in the left dot. Black (blue) lines show the results without (with) SO coupling. The labels S and 
r+ correspond to eigenstates without SO coupling. 
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FIG. 11: Probability that both the electrons occupy the same 
dot with and without SO coupling in the ground state and 
first excited state. 
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FIG. 12: Occupation of the single-electron orbitals of definite 
spatial parity and spin for two-electron ground state with (b) 
and without (a) SO coupling for the electron pair in the double 
dot. 
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FIG. 13: Contributions of the two-electron orbitals to the 
ground state with (b) and without (a) SO coupling for the 
electron pair in the double dot. 



ties for the discussed transitions - which are all forbidden 
in the absence of SO coupling - remain very small (less 
than 0.5%) as compared to the ones found for the single 
and three electrons. 

For two electrons the role of the electric field for the 
low-energy optical absorption is different from the single- 
electron case. For N = 1 the electric field distinctly shifts 
the energy of the absorption lines (Fig. 6). For N = 2 
the energy shift is very weak, only the transition prob- 
abilities are affected. For the single electron the energy 
shifts resulted from spatial electron-charge displacements 
of the initial and final states induced by the electric field. 
For two electrons these shifts are hampered (see Fig. 1) 
since the charge shift implies appearance of a double oc- 
cupation of one of the dots. Fig. |9] shows the charge 
localized in the left dot in function of the electric field. 
For 1 (and TV = 3) the dependence of the charge on 
Fx is the strongest at zero electric field, while for N = 2 
we find a plateau centered at = 0. 

For i? = we did not find any SO coupling infiuence 
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on the charge distribution as a function of the in-plane 
electric field. Nevertheless, such an effect is observed in 
the presence of the external magnetic field - see Fig. [TOl 
For B = 0.4 T the ground-state without SO coupling 
corresponds already to the spin triplet, in which - due 
to the Pauli exclusion - localization of both electrons in 
the same dot requires occupation of an excited single-dot 
energy level. The charge of the two-electron system for 
the triplet ground-state is even more resistant to shifts by 
the electric field than for the singlet state [compare Fig. 
fTHTb) and (d)]. For S = 0.4 T the ground-state becomes 
singlet again near 0.4 kV/cm. The electrons in the singlet 
state occupy more easilj^ the dot made deeper by the 
electric field which restores the singlet ground-state when 
Fx is applied. We notice [see the dashed line in Fig. 
[TUf d)] a jump in the occupation of the left dot at the 
singlet-triplet transition. For B = 0.6 T a similar effect is 
observed only at higher [the dashed line in Fig. fTUTf)]. 
The SO coupling mixes the singlet and triplet states and 
we notice that the electron charge in the left dot [blue 
lines in Fig. [TUTb.d.f)] becomes a smooth function of 
Fx- As a general rule, when the ground-state without 
SO coupling is singlet (triplet) - the SO coupling reduces 
(enhances) the occupation of the deeper dot. 

At the singlet-triplet transition the SO coupling in- 
fluences also the probability of finding both the elec- 
trons in the same dot (Fig. [TT|) . Without SO coupling 
the ground-state probability exhibits a rapid drop at the 
singlet-triplet transition near 0.4 T. The spin-orbit cou- 
pling infiucnccs the double occupation probability only 
for non-zero B. 

In order to quantify the occupation of the single- 
electron even and odd parity orbitals we first project 
the two-electron eigenstates of operator (11) into the ba- 
sis composed of single-electron eigenfunctions obtained 
without SO coupling (denoted as ip' in the following). 
For a state v we consider the projection in form 



dli = ^ E C'.';-<V'.(1)^,(2)-V.(2)^,(1)| 

i,j>i 

iv;(i)V';(2)-v^;(2)^;(i)). 



(13) 



An cigenfunction ip'f. has a definite spatial parity and 
z-component of the spin associated with a spinor Xk 
which is the Sz cigenfunction of eigenvalue fi/2 or ~h/2 
iXk = I t) or Xfc = I 4-))- Hence, the occupation of the 
spin-up even-parity single-electron wave functions can be 
calculated as 



=(et)= \dki\^Mk,+)Ss{kA) + Sp(l,+)Ss{lA)], 

(14) 



k,l>k 



where 



and 



Mfc,±) = l±IM(!rVi(z!0^ (15) 



The occupation of the spin-up odd-parity single-electron 
states is determined by the formula 

oc(ot)= E \dki?[5p{k,^)5s{kA)+5p{l,-)5s{lA)], 
k.lyk 

(17) 

with an obvious generalization for the spin-down compo- 
nents. The results are displayed in Fig. [121 Without 
SO coupling i) below 0.4 T the ground-state is even par- 
ity singlet - the electrons occupy mostly the even par- 
ity states ii) above 0.4 T the ground-state is odd parity 
triplet - the spin-down contributions are removed, one 
of the electrons occupy an even parity and the other an 
odd parity orbital. The jump of the occupations near 
0.4 T that is observed in the results without SO coupling 
is replaced by a smooth transition when SO coupling is 
applied. The values obtained for orbital occupations in 
both large and zero B limits are similar. 

Non-conservation of the spatial parity in the presence 
of SO coupling for the two-electron states becomes evi- 
dent when one considers contributions of the two-electron 
basis elements. The contributions of the elements in 
which both electrons occupy orbitals of the same spatial 
parity are calculated as 



\dki\^6p{k,+)6p{l,- 



k,l>k 



for the even parity orbitals and 

Coo = Y \dki\'^Sp{k,~)Sp{l,-), 



(18) 



(19) 



k,l>k 



SsikA)^{xk\t)- 



(16) 



for the odd parity orbitals. Contribution of the two- 
electron basis elements in which the electrons occupy op- 
posite parities is 

Coe+eo = Y \d-ki\'^ [Sp{k,-)Sp{l,+) + 5p{k,+)Sp{l,~)]. 
k.l>k 

(20) 

The results are displayed in Fig. [131 Without SO cou- 
pling for B < 0.4 T the contribution of the basis elements 
in which the electrons occupy opposite parity eigenstates 
is zero. In the triplet ground state for B > 0.4 T the 
electrons are bound to occupy orbitals of opposite pari- 
ties. When the SO is present for B = there is a nearly 
10% contribution of basis elements in which the electrons 
occupy orbitals of opposite parities. The Coe+eo grows 
with the magnetic field, but it stays below 100% in the 
studied range of B. This result and the ones presented 
above indicate that for two electrons the SO coupling has 
a noticeable infiuence on the ground-state properties in 
contrast to the single electron case. 



C. Three electrons 

For = 3 in the absence of SO coupling the magnetic 
field leads to the ground-state spin-polarization transi- 
tion near 3 T in both the single (Fig. [H)) and double 
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[Fig. [T5ja)] dots. For symmetric dots this transition is 
associated with energy level crossing even when SO cou- 
pling is introduced since the ground-states on both sides 
of the transition correspond to opposite s-parities. The 
in-plane electric field opens an avoided crossing at the 
ground-state spin polarization transition [see Fig. I17f a)]. 

For three electrons in a single dot without SO coupling 
one observes (Fig. [T4)) crossings of three s-odd energy 
levels near 2 T. For the double dot [Fig. [TSl a)] the cross- 
ings appear in more separated magnetic fields. The three 
crossing levels have different z projections of the spin. 
Similarly as for A'^ = 1 the SO coupling opens avoided 
crossing in the absorption spectrum, but for = 3 three 
energy levels participate in this avoided crossing instead 
of two. These avoided crossings arc accompanied by a 
smooth variation of the spin [Fig. [TSTb)]. 

For the spin unpolarized ground-state (B < 3 T) the 
lowest-energy optical transition goes from the ground- 
state to the odd s-parity states with (s^) ~ Tijl. Without 
SO coupling and in terms of occupation of single-electron 
orbitals we observe (Fig. [T6|) a transition of one of the 
electrons occupying a bonding orbital to an occupied an- 
tibonding orbital. One finds a single bright line similar 
to the one found for TV = 1. For i? > 3 T the principle 
line in the ground-state absorption spectrum disappears 
due to the ground-state spin polarization. 

The in-plane electric field increases the energy splitting 
between the ground-state and the first excited state lead- 
ing to an increase of the energy absorbed at the optical 
transition [Fig. [mb)]. The form of the avoided crossing 
is not affected by the field - like in the single electron 
case. 
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FIG. 14: Three-electron energy spectrum for the single elon- 
gated dot. Black (blue) lines show the results without (with) 
SO coupling. Solid (dashed) lines correspond to even (odd) 
s-parity states. The arrows in the plot indicate the z compo- 
nent of the spin without SO coupling. 



For the lowest-energy even s-parity state both occu- 
pation of single-electron spin-orbitals [Fig. [TSla.c)] and 
contribution of three-electron basis elements of definite 
spatial parity [Fig. [TW a.c)] are only weakly affected by 
both the magnetic field and the spin-orbit coupling. The 




FIG. 15: Three-electron energy spectrum (a) for the double 
dot. Black (color) lines show the results without (with) SO 
coupling. Solid (dashed) lines correspond to even (odd) s- 
parity. The short arrows in the plot indicate the z component 
of the spin without SO coupling and the longer ones show 
the allowed optical transitions from the ground-state, (b) z 
component of the spin for the lowest even s-parity and three 
lowest odd s-parity energy levels. Type and color of curves 
for these states is adopted of panel (a). 



(b) V,= 10meV 




FIG. 16: Optical absorption spectrum for the three-electron 
system in the single dot (a) and in the double dot (b). Black 
(blue) dots correspond to SO coupling absent (present). Area 
of the dot is proportional to the absorption probability. 



dependence of the studied quantities on the magnetic 
field is more spectacular for the lowest energy s-odd state 
[Figs. WLh,d) and[Ili;b,d)]. Without SO coupling the 
lowest-energy s-odd level corresponds to even parity only 
between 1.9 and 2.8 T, hence the vanishing contribution 
of the even-parity three-electron basis elements outside 
this B interval. In the presence of SO coupling the con- 
tribution of the even-parity basis elements extends over 
the entire studied range of the magnetic field. 
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FIG. 17: (a) Three electron energy spectrum for SO-coupled 
double dot at F^, = (black dotted lines) and for Fx — 0.5 
kV/cm (blue solid curves). The arrow indicates the ground- 
state avoided crossing which is opened in presence of non- 
zero Fx- (b) Optical absorption spectrum for the SO-coupled 
double dot. Black (blue) dots correspond to Fx = {Fx — 0.5 
kV/cm). 
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FIG. 18: Occupation of the single-electron orbitals of definite 
spin orientation and spatial parity without (a,b) and with 
(c,d) SO coupling for the lowest energy s-even (a,c) and s- 
odd (b,d) state for three electrons in the double dot. 



IV. SUMMARY AND CONCLUSIONS 

We have presented a systematic exact diagonalization 
study of one, two and three-electron spin-orbit coupled 
systems in double quantum dots. We discussed the mix- 
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FIG. 19: Contributions of single-electron orbitals of a given 
symmetry to the lowest energy three electron s-even (a,c) and 
s-odd (b,d) states, with (c,d) and without (a,b) SO coupling. 



ing of the bonding and antibonding electron orbitals by 
the SO coupling. We investigated occupation of even 
and odd parity orbitals, the energy and optical absorp- 
tion spectra in crossed electric and magnetic fields as well 
as the electron distribution. 

For one and three electrons confined in a pair of iden- 
tical dots we found that the spin-orbit coupling only 
weakly affects the ground-state properties. A strong 
mixing of bonding and antibonding orbitals due to the 
spin-orbit coupling was found in the lowest-energy ex- 
cited states. 

In contrast to the odd electron numbers, for two elec- 
trons the spin-orbit interaction affects the properties of 
the ground-state since the spin-polarization becomes a 
smooth transition instead of an abrupt singlet-triplet 
transformation. On the contrary, the spin polarization 
of the three electron system in symmetric dots is not 
affected by the spin-orbit coupling since the low- and 
high-spin ground-states correspond to opposite s-parities. 
For three electrons the SO coupling makes the spin- 
polarization continuous only when the confinement po- 
tential contains an in-plane asymmetry, e.g. introduced 
by an electric field. 

For odd electron numbers the spin-orbit-coupling- 
induced mixing of spatial parities of the first excited state 
opens characteristic avoided crossings in the optical ab- 
sorption spectrum. An in-plane electric field shifts the 
initial and final states of the optical transition to oppo- 
site dots. In consequence it distinctly increases the en- 
ergy of the optical transition at an expense of a reduced 
width of the avoided crossing. 

The low-energy optical absorption for two electrons is 
only allowed by the SO coupling. For two electrons the 
in-plane electric field lifts the spatial parity selection rules 
but docs not essentially perturb the energy of the optical 
transitions. 
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